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ABSTRACT 


Starting with the basic fatigue life formula 



an Inversion formula 1s derived In form 
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where z is a function of strainrange and the ratio c/b. 

The inversion formula is valid over the entire life range of enaineering 
interest for all materials examined. Coi.formity between the two equations 
is extremely close, suitable for all engineering problems. 

The approach used to Invert the life relation Is also suitable for the 
inversion of other formulas involving the sum of two power-law terms. 


INTRODUCTION 

The availability of a formula expressing life directly in terns of 
strainrange and mean stress is a great convenience in cumulative damage 
analysis. Data for such an analysis usually appear as a sequence of loading 
each characterized by a strainrange and mean stress, and the analysis proceeds by 
establishing t^e life for each loading condition, and suitably summing cycle 
ratios based on the number of cycles actually applied at each loading 
divided by the life that would occur if that loading persisted until 
failure occurred. Unfortunately, however, the life relationship for each 
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loading condition Is formulated by an expression containing two terms with 
life raised to a fractional exponent. No exact c1osed>form relaticm for 
life can be obtained because the transcendental equation does not lend 
itself to exact solution. 

A manber of approximate methods have been obtained In the past for 
approximate, although accurate. Inversion of the life relationship 
to obtain a direct expression for life N^. Ttiese are discussed In 
Ref. 1. In one form the equation, defined in Ref. (1) as the Collocation 
method. Is recast to an entirely different form, although numerically quite 
close to the actual life relation to be Inverted. However, the manerical 
conformity can be achieved over only a few decades of fatigue life; therefore 
a "floating" relationship Is required, changing the range for each application 
to Insure that the life involved falls within it. While the procedure lends 
itself to easy computer programming, it would obviously be preferable to 
avoid the need for such a "floating" system. 

A second approach, described In Ref. 1 as the spline point method, 
takes a step toward simplification by providing the Inversion relationship 
in just two analytical expressions, one applicable below the transition 
strainrange (the strainrange where the elastic and plastic coir^onents are 
equal), the other above the transition strainrange. While quite accurate 
and simple to program, the Inconvenience of having to account for a two- 
part analytical expression Is still present In this method. The authors 
have therefore continued the search for a single-expression closed-form 
relationship suitable over the entire life range of Interest In common 
engineering problems. 


3 


In this report we draw on the experience gained In Ref. 1 with the 
development of the two-part Inversion expression to establish such a single 
expression. The report describes the basis of the method, and Its applicability 
to a large number of materials comonly used In engineering design 


METHOD 

The Life Relation to be Inverted 

The form of the basic life rel.-'l.onship was first proposed by Manson 
(Ref. 2). Later the same basic equation was expressed by Morrow (Ref. 3) with 
a new notation which Is now comonly used 

Y = 4 + T 


Here Ac = applied strainrange 
= cycles to failure 

and 0 ^ = material constants designated as ductility 
coefficient and strength coefficient, respectively, 
b and c * material constants designated as the elastic 
and plastic exponents, respectively. 

An alternate form of the life relation has been expressed by Manson (Ref. 4) 


R 


t 



( 2 ) 


where R 

e 



Ac ■ applied strainrange 
■ cycles to failure 
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where and acj are transition life and strainrange given by 

:/(c-b) 




(3) 


«T ■ J 


l/(b-c) 


W 


Modification of the relation to account for mean stress Is discussed In 
[5]. In Ref. [6] It Is shown that, basically, the form of the relation Is 
still given by Eq. (2), except that: 

1. Nj is replaced by Nj, where 



where « mean stress 

2. The transition strainrange is the sane as for completely 
reversed loading, but can be replaced by k' if experimental Information 
Is available to Indicate that the cyclic stress-strain curve Is affected 

by mean stress. In the absence of such explicit Information Is taken 
as unity. 

3. The term In Eq. (5) may be replaced by k^j^o^ If experimental 
evidence exists to provide a more accurate relation between cyclic life and 
mean stress [7] than that originally proposed by Morrow in [8] that mean 
stress can be accounted for by replacing In Eq. (1) by - o^. 

Thus, whether mean stress Is or Is not present, the same basic equation 
of the form of Eq. (2) applies. We shall therefore focus our attention on 
the Inversion of this relation, recognizing that N.^. definitely depends on 
mean stress and ma^ be a weak function of mean stress. 
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Choice of the Form of the Inverted Relationship 

An Intuitive choice of form of the Inverted equation Is 



1 1 

R ^ + R ^ 
c c 


( 6 ) 


cq. (6) has the sane asymptotes as Eq. (2)« so ^ ut It Is a valid 
representation of Eq. (2) at extremes of life and strainrange, but 
unfortunately the two equations do not coincide for Intermediate values 
In the vicinity of the transition point for the parameters associated with 
most materials. Fig. 1 shows, for example, the results for T1-6A1-4V. The 
continuous curve Is the proper life relation obtained by adding vertical 
strain ordinates at any value of 2N^, while the dotted curve Is determined 
from Eq. (6), which essentially adds horizontal abscisas at any value of 
strainrange. Of course the simplified Inversion formula is correct at life 
extremes where the coordinate values on one line are negligible compare'' to 
those of the other; In the region of the transition point, where the 
coordinate values associated with both elastic and plastic lines are 
numerically significant relative to each other, a considerable difference 
exists between the results of Eq. (2) and those of Eq. (6). 

The next logical step Is to Introduce the parameter z according to 

Again It Is clear that the asymptotes of Eq. (7) are the same as those of 
Eq. (2), since they are 


1 1 
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HOMever. we now have an adjustable parameter z which can be chosen to 
produce good fit over a large range about the transition point while retaining 
the desired asymptotes for all choices of z (except the trivial case of z * 0). 
D atemlnatlon of the Parameter z 

The first approach was to determine whether z could be taken as a single 
constant. In this case only a single point on the curve could be force- 
fitted exactly. Consider, for example, the choice for fit at point (2) in 
Fig. 2 immediately above the transition point. Here R » 2 (since elastic and 
plastic strainranges are both equal to acj), and ■ 1. For this point, 

Eq. (7) becomes 

1 » [2^^^ + 2^^^] (9) 

This equation cannot be solved exactly, but a very good approximation can be 
obtained using the experience gained in Ref. 1. For example, by re-writing 
the equation In the form 

z c i 

1 , [2 c ^ ^ 2 ^} (9a) 

It becomes clear that since the only two parameters In this equation are 
z/c and c/b, that the equation Is really a statement z/c » f(c/b). Thus, 
choosing a series of values of c/b and solving numerically for z/c from 
Eq. (9a), we can plot the resulting values of z/c vs. the chosen c/b. 

A linear relation results when the plot Is on log-log coordinates, leading 
to 

, -.632 

Z 2 • -1.117 c(|) (10) 
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Another logictl point to consider for exact fit is point (1), where ■ 1, 

N A 

The value of ^ at this point has already been determined in [1] as e 

"t 

«t»re 1^1 1 (”) 

Substituting these values into cq.' (7) and solving by the same technique 
as used for solving Eq. (9) results in 




-0.889c 


It I 


>0.36 


02 ) 


Finally, a third point (3) in Fig. 2 can be used for fit. For this 
point is as much lower than as it is higher than N.j. at point ( 1 ). 
In other words, N, is the geometric mean of the life values at 
(3) and (1). Thus, since |ij- ' e at this point, use of Eq. (2) results 
in » e" + e"*' . Again substituting these coordinates into Eq. (7) 
and solving numerically, results in 


23 = -).Z37c(|) 


r.832 


(13) 


Any one of the three values of z from Eos. (10), (12), or (13) would provide 
a first-approximation constant value of z that can be used in Eq. (7). 

Fig. 3 shows the types of fit that can be achieved by using successively 
points (1), (2), and (3) for the particular alloy Ti-6A1-4V. While all 
fits at'e quite good, they are not adequate for applications reouiring high 
accuracy over the entire life range. By taking z as a second-degree 
polynomial in InR^, however, rather than a single constant, three adjustable 
constants become available, permitting fit at all three points ( 1 ), 

(2), and (3), and producing, therefore, an exceptionally good fit over the 
entire life range. The choice of location of the three points is exceptionally 
fortuitous, since it permits the final result to be expressed literally 
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In general terms of all the narameters Involved. Details of the procedure 
are outlined In Appendix A, the final result being 


where 


exp[P(ln RJ^ ♦ 0(ln R^.) +S] 

(14) 

c 2 c 

-.001 277 (|) + .03893 (^) - .0927 

(15) 

c 2 r 

.004176 (|) - 0.135 (|) + 0.2309 

(16) 

In J-.889C (|) 

(17) 


Final Formulation . The basic procedure Is, therefore, as follows: 

Starting with the commonly used form of the life relation Eq. (1) 
calculate the transition point coordinates from Eqs. (3) and (4), casting 
the basic relation according to Eq. (2). The Inverted relation Is Eq. (7) 
where z is given by Eq. (14) with the values of P,Q, and S given In Eqs. (15)- 
(17). If mean stress Is present, use the same procedure, except use 
Eq. (5) to obtain thr transition life. While In most cases the mean stress 
multiplier and transition life multiplier k^ can be taken as unity, actual 
values of these multipliers can be used If experimentally determined. All the 

equations necessary for the inversion procedure are summarized In Table II. 

DISCUSSION 


While the inversion procedure described in Ref. 1 Is more compact 
than that described In this report. Its drawback Is that It requires 
two formulas, one valid above the transition strainrange, the other 
below. The inversion procedure described here Is a single-expression 
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relation valid over the entire. life and strainrange of the material. 

Several smnewhat lengthy formulas are Involved In the procedure, but 
mmerlcally there Is little difficulty Involved since they can easily be 
programmed, even with simple hand>held calculators. The degree of 
conformity between the basic life relation and the Inverted relation Is 
remarkably good, conforming to any reasonable engineering standard likely 
to be required. 

Table I has been prepared to determine the conformity for a large 

nimiber of materials over the practical life range frwi 10 to 10® cycles to 

failure. The results are not shown for lives lower than 10 cycles or greater 

than 10® cycles, since agreement becomes even closer as the points Involved 

lie closer to the asymptotic elastic and plastic lines. Data for the table 

were obtained from Landgraf et al's compilation [9]. For each material the 

basic fatigue parameters are listed In columns 2 to 6, and the calculated 

transition point coordinates are shown in columns 7 and C. The strainrange 

required to produce lives of 10 to 10® cycles. In decade Increments, were 

first calculated from the basic life equations of form as In either Eq. (1) 

or (2). These strainranges were then used In the inversion formula, 

Eq. (7) to obtain the lives that the formula would yield, as c(^ipared to the 

2 

"exact" values 10, 10 , etc. These lives are shown In Colimms 9 to 14. 

An entry of unity means that the Inversion formula gives exactly the same 
results as the basic equation. The degree of discrepancy between the 
"exact" equation and the inversion formula is indicated by the departure from 
unity of the number entered at each life level. In most cases the error 
Is only two or three percentage points In life , a degree of accuracy far 
exceeding the experimental scatter usually associated with the detenninatlon 
of the basic life relation. Only one entry In the entire table Involves a 
difference of more than lot - material no. 1 at the 10 cycles life level. In 
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most cases the conformity Is extremely close. 

Figure 4 shows the life relations according to Eqs. (2) and (7) for the 
four nmterlals for which the discrepancies were the greatest among those listed 
In Table I. Even for these materials the difference between the two curves can 
not be discerned on a reasonable graph scale, except for the very lowest lives. 
Still, the discrepancies are negligible for engineering purposes. It can reason 
ably be concluded that the Inversions formula Is sufficiently accurate for all 
materials for all engineering purposes. 

The monotonic and cyclic stress-strain curves represent other examples of 
relations that can be Inverted by similar procedures. Appendix B outlines the 
application to the cyclic stress-strain curve. Since st**ess is expressed 
directly in terms of strain, closed form analytical expressions can easily be 
obtained for stress in terms of strain for double amplitude stress strain curves 
for both increasing and decreasing directions. With the proper rule for 
incorporating memory, it is possible to find the str jS response of a material 
for a given. complex strain histor'y. 

The basic inversion approach has utility for other applications, as well, 
whenever one variable is expressed as the sum of two negative power-law 
expressions of a second variable. The Smith-Watson-Topper [10] relationship for 
treating mean stress effects also involves such a formulation, and is therefore 
amenable to treatment by the method of this report. Appendix C presents the 
basic outline of the approach. Other potential applications are creep strain 
analysis, wherein a creep or creep-rupture curve is sometimes represented by 
the sum of two power-law relations. Numerous other applications can be 


envisioned. 



n 


CONCLUDING REMARKS 

The method presented in this report for inverting the life relationship 
provides extremely accurate results for all materials examined over the entire 
range of fatigue lives of interest in common engineering applications. The 
formulas involved can easily be programmed on a computer or hand-held calculator. 
No trial -and-error calculations are involved. All the relations involved are 
summarized for convenience in Table II. The basic inversion method may also 
find utility in other applications such as inverting the stress-strain curve, 
analytical treatment of alternative methods of calculating mean-stress effects 
on fatigue life — such as the Smith-Watson-Topper parameter, and inverting 
creep and creep-rupture relations involving two power-laws. The approach is 
even suggestive of how to treat relationships involving three or more power-law 
terms, although details would require further study. 


TABLE I 

APPLICATION OF INVERSION FORMULA FOR FIFTY CHARACTERIZED MATERIALS OF ENGINEERING INTEREST 
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TABLE II 


LIST OF FORMULAS USED IN INVERSION METHOD 


For completely reversed cycling: 
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^ (2N^)^ + c^(2N^)^ 

Transition strain-range Acj « 2(e^) * 


c/(c-b) 


"f 

Tj 


Transition life Nj = 


If the mean stress of is present, transition life 


I _ 1 r / OKI \"b ^0 -i-l/b 
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Inverted relation: 


"f ' "t • 


Z Z 
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where 


z = exp[P In^R^ + Q InR^ y S] 


P = -.001 277 (|) + . 03693 (|) - .0927 

Q = .004176(|)^ - .135(|) + .2309 

f c 

S ' lnl-.889c(|) 





Fig. 3. Exact and inverted life relations 
for Ti-6A1-4V {170®F). 


■ exact Equation 

-Fit throuqh point (1), Eq. (iZ) 
-Fit throuqh point (2), Eq. (10) 
.Fit throuqh point (3), Eq. (13) 


(I) 


I(A 

- El 

Fig, 4. Comparison of actual and inverted life relations for above four materials which 
show maximum discrepancy as shown in Table I. 
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Appendix A 


DETERMINATION OF THE Z-PARAMETER 


To fit z through three known points represented by Eqs. (10) » (12), 
and (13), we start by noting that the values of z/c for each of these equations 
Is a function of only (|). Furthermore, the values of R^ at points (1) 
and (2) are constants, while at point (3) It also depends on only c/b 

(sincefi- (|)°‘^ fron Ref. [1], e‘^® » exp^o.78 (|)“°‘®^1 and 

-cfi r c 0.3^ -* 

e ■ exoh].78 {■£) J, thus both terms are functions of only c/b). As a 

result. If we formulate the relation 


In I * P (In R^)^ + Q(1n R^) + s' 


(A-1) 


It Is clear that satisfying the three equations at point (1), (2) and (3) 
of Fig. 2 will produce relations for P, Q and S that will be functions 
of only c/b. However, since these functions are rather con^llcated because 
of the exponentials involved, we proceeded in more simplified fashion. 
Choosing specific values of c, c/b for actual materials, we applied the 
coordinates expressed In Eqs. (10), (12) and (13), to substitute Into 

s I 

Eq. (A-1) and solved for P, Q, and S . The value of 5 could easily be 

I 

determined In closed form, since at ■ 1, In R^ » 0, so S depends 
only on c/b, resulting in 


S 



(A-2) 


However, for determination of P and Q the graphs shown in Fig. A-1 were 
used. Calculations for P and Q are shown for various materials, and are 
plotted against (c/b). As expected, a single curve results. Passing a 



second degree polynomial through each of these two curves, using least-squares 
analysis resulted In the relations shown In Eqs. (15) and (16). By combining 
*he 1n^ term of Eq. (A-1) with the S' tern of Eq. (A-2), the resulting 
value of S shown In Eq. (17) was obtained, leading to the equation for z In 
Eq. (14) as a direct consequence of Eq. (A-1). 

Although the above method was adopted for final use. other approaches 
were Investigated during the study. Among then was one In which expressions 
for z similar to Eqs. (10). (12) and (13) were written for a nun^r of 
additional points along the curv>*. Least squares analysis was then applied 
to obtain a fit of an equation similar to Eq. (14) through the redundancy 
of points. The resulting equations were less accurate than the 
method finally adopted. Other approaches were more cwipllcated but not 
more accurate. The method adopted provides a high degree of accuracy with 
a relative simplicity of underlying basis and ease of final application, 
as discussed In the report. 



Fig. A-1. Determination of Constants P and Q in Eq. (A-1) for 
various materials shown. 
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APPENDIX B 


INVERSION OF Tff CYCLIC STRESS-STRAIN CURVE 
Tho cyclic stress-strain curve can be expressed by the following: 



(B-1) 

(B-2) 

(B-3) 


Nf 

iir-) can be expressed In terms of R, by Eq. (7). 
^ b 


Lo • EiiCj 




Here Eq. (B-3) becomes, 
(B-4) 


where 2 Is computed by Eqs. (14) to (17). 

Fig. (B-1) shows the basic cyclic stress-strain curve for T1-6A1-4V, and 
the Inverted relation by Eq. (B-4), which are essentially Identical. Hence 
stress Is known directly In terms of applied strain-range. 




Fig. B-1 


Comparison of basic cyclic stress-strain curve and 
inverted relation for Ti-6A1-4V (170®F). 


APPENDIX C 


INVERSION OF SMITH-WATSON-TOPPER RELATIONSHIP 


Life associated with a given hysteresis loop with mean stress can be computed 
by the knowledge of maximum stress and strain range by the SWT procedure [7]. 
According to this method a universalized curve results when we plot vs. 

for all combinations of Ae and mean stress. Thus If we start with case for 
conpletely reversed loading. 


0 

max 



Using Eq. (1 ), we get 



°max^ 


0 ^ 


(C-1) 


Expressing Eq. (C-1) in the form of Eq. (2), 
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where 

"t 

Is given by 

Eq. (4) 






m m 

2b 
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Eq. (C-2) can be easily Inverted by using the procedure described for 
Inverting strain-life relationship 



2 


P * 

Q = 
s = 

Actual (f E 


as shown in Fig. 



loop, if o_,^ and 

iTIqX 


exp[P In^ R + Q In R ♦ S] 

2 

-.001277[^] + .03893[^] -.0927 
2 

.004176 [^] - .135[^] + .2309 
ln|^-.889(c+b)(^) 

) vs. curve and inverted curve for Ti-6A1-4V 
!)-l are essentially indistinguishable. For a given material 
P,Q,S can be easily calculated. Then for any given hysteresi 
^ are known, can be determined by Eq. (C-4). 






